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1. Introduction 

The rare decays of B mesons have never been of greater interest, both experimentally 
and theoretically. The first observation of a decay mediated by the quark transition b ^ s 
recently has been reported by the CLEO Collaboration [Q , who found a branching fraction 
for the process B K*'y of (4.5 ± 1.9 ±0.9) x 10~^. Such transitions are typically induced 
by the exchange of virtual heavy quanta, the effects of which appear at low energies as local 
operators multiplied by small coefficients. It is hoped that the detection of these suppessed 
interactions in the guise of rare B decays may provide a direct window to physics at much 
higher scales. 

In order for such a hope to be realised, however, it is necessary to connect the quark- 
level operators which are generated perturbatively to the hadronic transitions which are 
actually observed. This involves the consideration of nonperturbative hadronic matrix 
elements, which typically are incalculable. One common approach to this problem is to 
consider inclusive rates such as B ^ Xg rather than individual exclusive channels, and 
to model the inclusive transition by the decay of a free bottom quark to a free strange 
quark. It is hoped that for the b quark mass mi, sufficiently large, the operator mediating 
b — i> s acts over distances short compared to the scales of confinement and strong QCD 
interactions, and the approximation is a good one. 

The issue of how good this approximation really is originally was addressed by Chay, 
Georgi and Grinstein 0. Using the tools of the heavy quark effective theory (HQET), 
they showed that the free quark model is in fact the first term in a controlled expansion in 
l/m^, and hence is arbitrarily accurate as rrib ^ oo. In addition, they demonstrated that 
there are no contributions to the rate at subleading (1/mb) order, and that any corrections 
could only come in at order or higher. 

In this paper, we extend the work of Chay et al. to compute the leading corrections 
to free quark decay, for the inclusive processes B — ^ and B — > Xs(-^l~ . While 

the corrections here are not particularly large, it is important to know their size if 

the free quark decay model is to be trusted. We also believe that our computation is a 
very nontrivial application of HQET in a somewhat unfamiliar regime, and is hence quite 
interesting in its own right. 

Finally, we note that work which overlaps with ours has been performed recently, in 
a somewhat different formalism, by Bigi et al. 0]. 
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2. The Operator Product Expansion and Matrix Elements in HQET 

In this section we will discuss our procedure in general terms, to elucidate the structure 
of the expansion before beseiging the reader with particular details. We are interested in 
the rare decays of b quarks, such as b — S7 or b — > se+e~, which are mediated at low 
energies by local operators of the form 

OP{(l)) = sThP{(p) . (2.1) 

Here P{(p) is meant to stand for some function of perturbatively interacting fields such 
as leptons or a photon, and F is a general Dirac structure. Interactions such as ( |2.1|) are 
typically induced at high energies by the exchange of virtual W bosons, top quarks, or new 
exotic quanta. At low energies they appear in the effective Hamiltonian as local operators, 
with coefficients which may be computed using renormalisation group techniques. We will 
take the presence of such operators simply as given; our interest will be in the evaluation of 
their hadronic matrix elements. We note that operators of the form ( |2.1| ) are not the only 
relevant ones which will appear at low energies; for example, we will typically find four- 
quark operators as well. For these, the techniques which we will present below will only be 
appropriate when the invariant mass of the intermediate qq pair is far from any quarkonium 
resonances. We will return to this issue in our discussion of the decay b — se"'"e~. 

For now, however, we restrict ouselves to operators with the structure ( |2.1| ). They 
induce quark level transitions of the form b — s. However, since the quarks are confined, 
what is observed is the decay B Xg , in which a B meson decays to an arbitrary hadronic 
state Xs with strangeness 5" = — 1. (Decays from the lowest lying bottom baryon, Af,, are 
also possible.) Hence we need to compute matrix elements of the form 

{X,---\OP{c/>)\B), (2.2) 

where the ellipses denote denote the additional perturbatively interacting fields which cou- 
ple to P{(p)- Unfortunately, exclusive hadronic matrix elements such as ( p.2| ) are governed 
by nonperturbative strong interactions and are typically incalculable. At best, SU(3) 
and heavy quark symmetries may be used to relate the form factors which appear in one 
such matrix element to those which appear in another But computations from first 
principles are not at this point possible. 

Instead of considering the exclusive modes individually, then, we will we will sum over 
all possible strange final states Xs- As has been shown by Chay et al. [||, the inclusive decay 
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rate may in fact be calculated reliably. Previous computations of the inclusive rate have 
relied on the free quark decay model, in which the sum over exclusive decays is modeled by 
the decay of an on-shell bottom quark to an on-shell strange quark. For nib oo this is 
justified by arguing that the decay is essentially a short distance process, which occurs on 
time scales much shorter than those which govern the eventual hadronisation of the final 
state. This argument can be made precise within a controlled expansion in inverse powers 
of the bottom quark mass 0, and we will be able to compute the leading corrections 
to this limit. 

Squaring the matrix element ( |2.2| ) and summing over Xg, we find a differential decay 
rate of the form 

= ^ $^d[P.S.](27r)45W(P5 - - q){B\ iOtp(</,)t \X, ■ ■ ■) (X, ■ ■ • | iOP(</>) \B) . 

(2.3) 

Here Pg and Px are the momenta of the intial B and final Xs systems, and q = Pb — Px 
is the momentum transfered to the other decay products. The symbol d[P.S.] denotes an 
appropriate phase space differential. The part of dF which involves the fields P{(p) may be 
calculated perturbatively. We then find that dF is equal to the product of known factors 
times an expression W{q) which involves only the quark and gluon fields: 

W{q) = J](27r)4(5(4)(PB -Px- q){B\ O \B) . (2.4) 

Xs 

Here the sum over Xg includes the hadronic phase space integral. The treatment of this 
nonperturbative expression is the subject of the rest of this section. 

We begin by noting that W{q), being essentially a total decay rate, is related by the 
optical theorem to the discontinuity in a forward scattering amplitude. That is, we may 
write 

W{q) = 21mT{q) , (2.5) 
where an example of the time-ordered product 

T{q) = {B\T{OlO}\B) (2.6) 

is shown in fig. |I|. 

Now we come to a crucial observation [Q. The sum over Xg in eq. (^]^) includes 
hadronic states with a large range of invariant masses, < P^ < M'^. The energy 
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which flows into the hadronic system Xg scales with mi, as the bottom mass increases, 
and in the hmit rrib ^ oo is typicaUy much larger than the energy scale Aqcd which 
characterizes the strong interactions. Hence, in all but a corner of the Dalitz plot, in 
which P|- ^ m^, the strange quark in fig. |l| is far from its mass shell. In position space, 
this means that the points at which O and act must be very near each other on the scale 
of nonperturbative QCD, and it is appropriate to perform an operator product expansion of 
the time-ordered product in eq. (pT^). This operator product expansion may be computed 
perturbatively in as{mi,). It will be valid over almost all of the Dalitz plot, failing only in 
the region where P|- is small. In the large limit, the fractional contribution of this bad 
region to the total phase space integral is negligible, and our calculation of the inclusive 
decay rate based on this expansion will be reliable. Our approach, then, will be to perform 
a systematic expansion in inverse powers of mt, of which the leading term will be the 
result in the mi, oo limit of the theory 0. However, we will also be able to compute 
the leading corrections to this limit, using the tools of the heavy quark effective theory. 

In this section we will discuss the form of the operator product expansion, and how 
to take the hadronic matrix elements of the operators which come out of it. When we 
apply this formalism in the following sections, the expressions which we derive sometimes 
will be quite lengthy. Here we will concentrate only on the structure of the procedure. In 
general, then, the time-ordered product ( |2.6| ) may be expanded in a series of local operators 
suppressed by powers of the mass of the bottom quark. 



T{(9t,(9}°^^J- 

mb 



Zmb 4m^ 



(2.7) 



The operator On is an operator of dimension 3 -|- n, with n derivatives. 

At this point, it is useful to introduce the heavy quark effective theory (HQET) 0], an 
effective theory of QCD in which the mass of the b quark is taken to infinity. This effective 
theory implements on the lagrangian level the new "spin-flavor" symmetry of QCD which 
arises in this limit [0. Both the mass and the spin of the b quark decouple from the soft 
bound state dynamics of the hadron of which it is a part; so far as the light degrees of 
freedom are concerned, the heavy quark is nothing but a static, point-like source of color. 
The exchange of soft gluons with the light degrees of freedom leave the b quark always 
almost on shell. Thus we can write its four-momentum p*^ as the sum of its "on-shell" 
momentum m^v^ and a "residual momentum" , such that the components of are 



always small compared to m^. It is then convenient to replace the usual quark field b(x) 
by a new two-component field h{x) with fixed four-velocity v^, 

h{x) = e''^''^-"P+b(x) , (2.8) 

where P-f = ^{1 + projects onto the quark, rather than antiquark, degrees of freedom. 
This effective field has the property that a derivative acting on h{x) yields the residual 
momentum k'^ rather than the full momentum p'^. An expansion in terms of D^/m;, then 
becomes sensible. Expanding in powers of l/rub, we may invert ( p.8|) to find 

b(x) = e"'"^""-^ [1 + ip/2mb + ■■■] h{x) . (2.9) 

Inserting this into the usual QCD lagrangian bipb, we find the effective lagrangian for 
HQET i, 

C = hv-im + 5C (2.10) 

where the correction terms [0 

5C= ^XiTifh- ^ZMKv-iTifh+ ^Z2{^i)hs^'''G^,h + 0{l/Aml) (2.11) 

are treated as perturbations to the rrib oo limit. Here the gluon field strength is defined 
by G^i, = [iD^,iD,^], and s^'^ = —^a^^. The renormalisation constants are given by 

/ / \ \ 8/25 

/ X N 9/25 ^ ' 



^2(^) 



above the charm threshold. 

Because the operator product expansion ( P77| ) is an expansion in D^/m^, we must 
express the operators On in terms of the HQET field h{x) rather than the full fields b(x). 
However, as we shall see, it turns out to be convenient to leave the leading operator in 
terms of b(x), and to expand the rest in h[x). The operators On which appear in the 
expansion (|2.7[) then take the form 

Co = brb, 

Oi=~hViD^h, (2.13) 
C2 = 7^riD^iD,/i, 
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and so forth. In each case, F denotes an arbitrary Dirac structure, in which we also absorb 
all dependence on the external momentum g, as well as on any other variables. We will 
keep operators in the expansion with up to two derivatives. 

We now turn to the evaluation of the forward matrix elements of the operators On 
between B meson states. At leading order, we need matrix elements of the form 

{B\hTh\B), (2.14) 

which is nonzero only for F = 1 or F = 7'^. In the second case, the conservation of the 
b-number current in QCD yields the matrix element normalised absolutely, 

{B\h'y^h\B) = 2P^. (2.15) 



This, of course, is why we left Oq in terms of the field b(x) in eq. ( p.l3| ). As for the scalar 



current, it may be rewritten in terms of the vector current plus higher dimension operators 
of the form of C2 1, 

bb = v;bj^h + [{iDf - {v ■ iD)2 + s^^G^^] h + .... (2.16) 



This identity may easily be proven by using eq. (p.9| ) to expand both sides in terms of the 
effective field h. It is only meaningful when the four-velocity f ^ of the b field is fixed. 
The correction term in eq. ( |2.16|) may be absorbed into 02- Hence, the leading term in 



the expansion of T{q) may be evaluated unambiguously, using eq. (|2.15|) . In fact, the 



leading term is precisely the free quark decay model result, which becomes exact in the 
limit nib ^ 00 0. The subleading operators On in the operator product expansion ( p7?| ) 
will provide systematically the corrections for finite b quark mass. 

The evaluation of the matrix elements of the higher dimension operators Oi and O2 
involves the equation of motion of the effective theory . This is given by the lowest order 
lagrangian, 

viDh = 0. (2.17) 
Since the external states are characterized only by their four- velocity v'^, Lorentz invariance 



severely restricts the forward matrix elements of operators of the form ( 2.13 ). For the 
operator Oi of dimension four, we find 

(M| Oi \M) = (M| hViBf^h \M) = (M| hVv^v ■ iBh \M) . (2.18) 
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However, this is now the matrix element of an operator which vanishes by the equation of 
motion ( |2.17| ). Pohtzer |^] has shown that aU such matrix elements vanish identically; his 
proof is outlined in the Appendix. Since Oi is the only possible source of corrections of 
order l/rrib to the lowest order result, we see that the leading corrections to the free quark 
decay model are actually of second order in the heavy quark expansion. As first pointed out 
by Chay et ai, this is a most surprising result, since exclusive decay modes all presumably 
receive corrections already at order l/rrib- Somehow these individual contributions must 
cancel in the inclusive rate. 

The dimension five operators do give nonvanishing contributions, of order How- 
ever, their forward matrix elements have a very simple parameterisation [jl^. The symme- 
tries of the effective theory may be used to write the matrix element as an ordinary Dirac 
trace, 

{M\02\N) = {M\hTiD^iD,h\M) = MBTT{TP+t(j^,P+} , (2.19) 

where 

V'/x:^ = ^>^li9^lu - v^v^) + ^A2io-^i. . (2.20) 

The mass parameters Ai and A2 are defined in terms of certain expectation values in the 
effective theory, 

(2.21) 

where (Im = 3 and (Im* = —1- The /U-dependence of A2 cancels that of the renormalisation 
constant Z2{p) ( p.l2| ). We note that Z2{mi,) = 1; hence from this point on we will drop it 
and by A2 mean A2(TOb)- 

The role which these parameters play in the effective theory is revealed when one 
expands the masses of the heavy pseudoscalar and vector mesons in powers of l/m^: 

Mb =mb + A- -^{Xi + 3A2) + . . . , 

(2.22) 

Mb* =mi, + A- - — (Ai - A2) + . . . . 
Zrub 

In this expansion, the term A represents the energy of the light degrees of freedom in 
the meson. We see that Ai and A2 are higher order effects of the finite b quark mass; 
Ai is essentially a "Fermi motion" effect, while A2, the leading spin symmetry-violating 



correction, arises from the hyperfine chromomagnetic interactions. From ( p.22|) we have 
the weU-known relation 



A2 = "^{Mb* - Mb) ~ ]{Ml* - Ml) = 0.12 GeV^ , (2.23) 
o 4 

where we are neglecting higher-order corrections in 1/mb- There have been attempts to 



extract both Ai and A2 from QCD sum rules by computing the matrix elements ( p. 21 
with the results [jll[ 

< Ai < 1 GeV^ , 

(2.24) 

A2(l GeV) = 0.12 ± 0.02 GeV^ . 
The parameter A2 is much better determined in this approach and agrees nicely with the 



experimental B B* mass splitting (|2.23|) . 

Finally, there is one other source of corrections of order l/m^, namely time-ordered 
products of Oi with the correction SC to the effective lagrangian ( |2.10|) . As discussed in the 
Appendix, these arise because the states \M) in the effective theory differ at order l/nii, 
from those \B) of QCD. The difference is compensated for at each order by computing 
matrix elements of the form |]Ty] 



i / dx{M\T{hrh,5C{x)}\M) , 

(2.25) 

dx {M\ T{hTiD^h , 5C{x)}\M) , 



and so on. We did not encounter the time-ordered product SC with the dimension three 
operator Oq above, because we were able to compute its matrix element ( |2.15| ) directly in 
full QCD. For the operators of dimension four, since we are working in the effective theory, 
we must evaluate eq. ( |2.25| ). First, we use the fact that the external states depend only on 
the four- velocity to write the analogue of eq. ( |2.18| ): 



i j dx{M\T{hTiD^h, 5C{x)\M) = i j dx{M\T{hTv^vim, 5C{x)]\M) . (2.26) 

We may now apply the identity derived in the Appendix, which exploits the fact that the 
operator which appears on the right-hand side of eq. ( |2.26 ) vanishes by the equation of 
motion of the effective theory. We then obtain the matrix element of a local current, 

i f dx{M\T{hTiD^h, 5C{x)]\M) = -v^- ^{M\hTP+F{T))h\M) , (2.27) 



where F(D) — (iD)^ + s'^'^G^^. This matrix element may then be evaluated in terms of Ai 
and A2, using eq. ( |2.19| ). 

The result of this long and involved procedure is an expression for the nonperturbative 
hadronic quantity T{q), of the form 

T(g) = To(g^ v q) + ^T2{q^ v q) + . . . . (2.28) 

The expansion is a series in l/mb and as(m;,), and the ellipses in ( p. 28 ) denote higher order 
terms in both small parameters. (The radiative corrections to To have been computed 
previously [l^]; we will not include those to T2.) We now take the imaginary part of T{q) 
to recover W{q), multiply by the perturbative part of the matrix element which couples 
to P{(p), and compute the inclusive differential width dF. This may then be integrated to 
give the total width, F, or other smoothly weighted distributions. As we have mentioned, 
the contribution of Tq will be precisely that of the free quark decay model . The leading 
corrections to the nii, 00 limit are of relative order and encoded in T2; they are 

expressible entirely in terms of the mass parameters Ai and A2. We will now apply this 
procedure, and compute these corrections, for two interesting examples. 



3. Application to Rare B Decays 

We will consider inclusive decays of the form B Xg'j and B Xsi~^i~ , where i = e 
or ^ is a light lepton. They are governed by the effective Hamiltonian density 

4Gf- 



Heff = --^VttVZj2^MOM ' (3-1) 



where the sum is over the truncated set of local operators 

01 = SaT^Piba C/37/,PlC/3 

02 = Sa7^PLb/J C/37^PLCa 

e 



1671"= (3.2) 
e2 _ 

Os = 7^ Sa7'^PLbae7/^e 



167r2 
167r^ 



O9 = 7^ Sa7''PLbae7/,75e . 



Here Pl = |(1 ~ 75) and Pr = |(1 + 75) are helicity projection operators, F^^, is the 
photon field strength and a, (3 are colour indices. We have included in Og and Og only the 
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coupling to the electron current; the coupling to the muon is analogous. There are also 
additional operators, such as Sa'^^PLha [u/37^Plu^ + . . . + b/37^PLb^], which contribute to 
these decays, but their coefficients are small and we shall neglect them. The coefficients 
Cj{mb) have been calculated in leading logarithmic approximation, both in the standard 
model and in certain minimal extensions, and are presented in refs. ||13|| -|[T9| 



We now apply the procedure of the previous section to compute the rates for inclusive 
decays mediated by the operators (|3.2[) . The first step is to construct the operator product 



expansion (|2.7|) , which takes the form 

1 „ 1 



T{bris,sr2b}°='' — 
nib 



(3.3) 



For now, we will allow Fi and F2 to be arbitrary Dirac matrices. To compute the terms 
in this series, we must expand the diagrams in fig. |l] in powers of l/nib- Fixing the four- 
velocity of the external b quark to be f ^, we may expand its momentum as = nibV^ -\-k^ . 
Then the graph in fig. ^(a) gives 

... ._ Fi (mf,^ -^ + ^ + ms) F2 

iM = -lUb ^ 2u^2 2^ 

[mbV — q + ky — + le 
UbVi ii) - $ + rhs) V2Ub 



^-^Ub ( -Fi^F2 - ^Fi - I + rhs)r2 {v - q)"koA Ub 

+ 0{l/ml) 

where q = q/rub, thg = Wg/m;,, and 

x = l — 2v-q + cf' — ml + ie (3-5) 

contains the pole corresponding to an on-shell strange quark, near the end of the physical 
cut. The spinor Ub which appears is the ordinary on-shell b quark spinor of QCD. From the 
matrix element (|3.4|) , we may deduce the first two terms in the operator product expansion 

my- 

00 = -bFi(^-| + m,)F2b, 

2- 4 - (3.6) 

01 = -/iFi7"F2iD„/i- — (t;-^)"/iFi(^-| + m,)F2iD„/i. 

Note that while we have left the leading operator in terms of the four-component fields 
b(x), we have expanded Oi in terms of the two-component effective fields h{x). The reason 
for this choice was discussed in the previous section. 
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To obtain O2 we must also expand the one-gluon graph in fig. |T](b), in order to identify 
the contribution from the gluon field strength G^i, = [iD^,iD,^]. Additional contributions 
to O2 arise when the full QCD fields b(a;) in Oiare replaced by the effective theory fields 
h{x) via the relation ( |2.9|) . Equivalently, one may expand the spinors ui, in the matrix 
element (|3.4|) in terms of the two-component spinors Uh of HQET, 



Ub 



1 + 



b 



1 



Uh 



(3.7) 



and check that the result may be made covariant. Finally, there will be corrections at order 
I /ml if the leading operator Oq contains a scalar current bb, because of the expansion 
( p.l6|) . These are still contained in Oq and are not included in 02- A straightforward 
calculation then yields 



X 



— hT,{f,-i + m,)T2{iBfh 

— {v - g)^7^ri7"r2 (iD^iD^ + iD^iD,)/i 



H — T^rhs hVi iaanT2G'^^h 



X 

2 



(3.8) 



x"^ 

+ (7'^ri7"r2 + ri7^r27") iD^iD^/i 

4 



x^ 



-qrhTi{i,-i + m,)r27^ iD^iD^j/i . 
To continue any further, we must specify the Dirac structures Fi and F2. 



3.1. B ^ Xs-f 

For the transition B —>■ Xg'j, only the operator O7 from ( |3.2| ) contributes and the 
operator product expansion simplifies considerably. We now contract the terms in the 
time ordered product (|3.3| ) with the external photon fields and take the matrix element 
between B mesons to construct the hadronic object T{q) defined in eq. ( p.6| ). Because the 
decay is to an on-shell photon, q'^ = is fixed, and T{q) becomes a function only of the 
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scaled photon energy v ■ q = E^/mb- Including the matrix element of the photon field, we 
find 

f{v ■ q) = i^{B\T{ha^''s,saP^h}\B) ■ {l{q, e)\ F^^F^^ |7(g,e)) 



-16MBmb{v ■ qf 



(3-9) 

1 Ai ( b — Qv ■ q\ ^ X2 3 



X 2ml V Sx-^ / 2ml -^^ 



The sum is over the transverse polarizations of the photon, and there is a factor of i 
from each insertion of the effective Hamiltonian (|3.1| ). We neglect contributions of order 
m^. Note that we distinguish between m^, the bottom quark mass which arises in the 
operator product expansion, and Mb, the B meson mass which arises from the relativistic 
normalisation of the states (and therefore drops out of the final expression). The inclusive 
rate for B Xg'j is then given by 



aGl ml 2i I m2t [ d^/c ~ 



(3.10) 

327r5 Mb 

where z = v ■ q, and the contour integral is taken around the pole at x = 0. It is 
straightforward to evaluate this integral, and we find 



(3.11) 



This expression for the total rate agrees with the result of ref. 0]. The first term is just 
what one would obtain in the free quark decay model. 

We may consider using the same method to compute certain features of the photon 
energy spectrum. Of course, the precise shape of the spectrum is not available to us, in 
particular its behaviour near the endpoint of maximum E^. This is true at any order in 
the 1 / mb expansion, because in this region the strange quark approaches its mass shell and 
the operator product expansion breaks down. For example, the fact that the true endpoint 
of the photon energy spectrum is found not at E^ = mb/2 but rather at (M^ — M]^) /2Mb 
is entirely unavailable to us in this formalism. 

It is instructive, however, to generalise eq. ( p.lO|) to calculate the nonperturbative 
contributions to the moments of the energy spectrum. For example, the deviation of the 
average photon energy from that of the free quark decay model is 



{E,) =Tr^^^\Vt,V:,W,{mb)\' ^ {mbz)zT{z) dz 
mb f ^ Ai + 3A2 
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2 V 2ml 



(3.12) 



This result has an interesting structure, if we compare it to the 1 / expansion of the B 
meson mass ( |2.22| ). Naively, we might have expected (E^) to be shifted from the free quark 
value mb/2 to half the physical meson mass Mb/2. However, that is not what we find; only 
the order 1/m^ terms contribute. The reason is that the correction A to Mb in eq. ( |2.22| ) 
is the contribution to the meson mass of the light antiquark and the other light degrees of 
freedom, which in this formalism are mere spectators to the decay. Since they are present 
in the final state Xg as well as in the initial state, they do not represent additional energy 
available to the photon. By contrast, the higher-order mass corrections proportional to 
Ai and A2 arise from terms in the effective Hamiltonian of the b quark, representing its 
"Fermi motion" and its chromomagnetic interaction with the soft hadronic surroundings. 
These bound state shifts in the b quark energy are then refiected in the average photon 
energy {E^). 

We could also generalise eq. to higher moments of the photon spectrum. How- 

ever, there arises an additional complication if we insist that the moments we compute be 
experimentally meaningful quantities. This is because they are constructed by convolv- 
ing a power of the photon energy with the measured energy spectrum, but this spectrum 
is only related to our computation once our result (|3.9| ) has been smeared over typical 
hadronic scales. That is, T{z) should be replaced by the smoothed quantity 



Tf{z)= j dz' f{z-z')T{z'), (3.13) 

where f{x) is some smearing function of width 5 = iS.E jmi,- If we take /(x) to be a Gaus- 
sian distribution, /(x) = exp(— V 7r5^, we can calculate the moments analytically. 
Keeping terms of order 5^, for the n'th moment we find 



{/\Ef n(n + 2) Ai ^ A2 



1 -f 2n(n — 1)— z — 2 7^ — 2 ~ ^""t 



2m^ 3 2m^ 2m^ 



(3.14) 



The total rate (n = 0) and average energy (n = 1) which we have already presented are 
unaffected by this procedure, but the same is not true for the moments with n > 2. Note 
that the effect of the smearing is proportional to 5^ rather than 5, and so is formally of 
the same order as the nonperturbative corrections we have been considering. However, 
in order for our inclusive predictions to be meaningful, the resolution Ai? with which the 
photon energy spectrum is measured actually must be much greater than Aqcd, so that 
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many exclusive states are always summed over. Hence, it is in fact this resolution, rather 
than the nonperturbative effects, which will dominate the corrections to the moments T^'^\ 
In addition, real gluon emission will broaden the energy spectrum over the entire allowed 
phase space < < mi,/ 2, which will affect substantially the shape of the experimentally 



measured spectrum [|T2 
3.2. B ^Xsi+i- 

The transition B Xsi~^i~ receives contributions from the complete set of operators 
in (|3.2| ). In particular, unlike the decay B ^sT, the four-quark operators Oi and O2 have 
non-vanishing matrix elements. In order for our treatment of the four-quark operators to 
be valid, it is crucial that the invariant mass of the lepton pair not be near any resonances 
in the charm system such as the 1/', so that strong final state interaction corrections will 
be small. In this case we can treat the contributions from Oi and O2 as effectively local 
on the scale of hadronic interactions. 

It is convenient to separate the total rate for B Xsi~^i~ into two terms, correspond- 
ing to the decay to left- and right-handed leptons 

dTB^x^e+i-^W^''iq^vq)Lf:, + W^"" {q\v q)L% , (3.15) 

where the lepton tensors are given by 

]^L{R)^,u ^ ^ pu^pt^_ _ gi^Up^ . p_ ^ ieA'^^Pp^^p_^ . (3.16) 

Here and p- are respectively the four-momenta of the and l~ . Since we are restrict- 
ing ourselves to £ = e or /x, we neglect the masses of the leptons. 

The two Lorentz structures which arise from the effective Hamiltonian ( p.2|) are 7^ ( 1 — 
75) and a ^1,(1 +75)1?!^. Hence it is convenient to write 



(3.17) 



j.L(R) _pL(i?)t 
-■-1 —^2 ' 

where s = g^. The A's and S's are then combinations of the coefficients ci, . . . , cg: 

= cs{mb) - cg{mb) + [3ci{mb) + C2{mb)] g{mc/mb, s) , 
= cs{mb) + cginib) + [3ci(mb) + C2{mb)] g{mc/mb, s) , (3.18) 
BL=b'^ = -2cj{mb) , 
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where the function g{mc/mi,, s) multiplying ci and C2 arises from taking the one-loop 
matrix elements of Oi and O2 and has the form 



= --log.^ + - + -^--^l-- 2 + - 



, V 1 - 4-2^/5 + 1 + ie 

X log = 

\v^l-4z2/s-l + ie^ 

Integrating over the lepton phase space, the total decay rate is given by 



(3.19) 



B^x^e+e~ i (27r)32£;+ (27r)32£;_ 



ml 



(3.20) 



We must next perform the contour integral in the v ■ q plane and then the Ei integral 
to obtain the differential decay width. Since the calculation is quite tedious and the 
intermediate expressions extremely lengthy, we present only the final result: 

dV j^^x .e+ e- G^^a^ 



ds 2567r5 

X 



- {1- s)Re{B*'A') 



i=L,R 



+ ^[{-¥' + ¥^l)\^'?~W + m^' (3.21) 



+ (s - f) Re (S*^A^)] 



+ ^[(-55^ + ¥5 + ^)1^1^ -1(1 + ^11^1^ 
+ (7s - 5)Re (S*M' 



The summation is over the two chirality states of the leptons. 

The leading term in eq. ( |3.21| ) reproduces the free quark decay model result obtained 



in refs. [jI^[|T8[, while the subsequent terms are the leading non-perturbative contributions 
to the decay rate. It is interesting to note that unlike the parton level result, which 
has a characteristic 1/s behaviour at small s from the one-photon intermediate state, the 
non-perturbative corrections approach a finite constant value as s ^ 0. The differential 
spectrum for the invariant mass of the lepton pair is plotted in fig. |2|. We have chosen 

16 



a top quark mass of rrit = 150 GeV, along with = 4.5 GeV, as{mw) = 0.12 and 
«s("^6) = 0.21, to generate the spectrum. The free quark decay model result (Ai = A2 = 0) 
is presented along with the spectrum for Ai = 0.5 GeV^ and A2 = 0.12 GeV'^. We have 
normalised the width for this decay to that for semileptonic B decay (which includes the 
nonperturbative corrections given in ref. [^). The modification to the B — ^ Xsi~^i~ rate 
is reasonably large and tends to enhance the overall rate for high mass lepton pairs by 
order 10%. 



4. Summary and Conclusions 

Because of the necessary cuts to remove backgrounds, the full spectrum from a decay 
such as S — > Xg'y and B Xsi~^i~ is not available in an accelerator experiment. It is 
therefore important to understand well the shapes of these spectra if one is to relate the 
observed branching fractions to fundamental parameters of the electroweak theory. This 
is particularly true for the high photon energy and high invariant lepton mass regions of 
the Dalitz plot. Modifications to the simplest model, that of free quark decay, arise from 
strong interactions that can be classified heuristically as perturbative and non-perturbative 
corrections. 

The perturbative corrections arising from gluon bremstrahlung and one-loop effects 
for B Xg'y have been computed previously fl^. It is to the non-perturbative corrrec- 
tions that we have addressed ourselves in this paper. We have detailed the formalism for 
treating the semi-hadronic inclusive decays of mesons containing a single heavy quark. 
Upon summing over all hadronic final states, one may express the rate for a given process 
in terms of a time-ordered product of quark bilinears. This time-ordered product is then 
expanded in a series of local operators, the matrix elements of which either are known 
or may be parameterised simply. Heavy quark symmetries and the heavy quark effective 
theory play a key role in the analysis. 

We have applied these tools to the rare decays B Xg'^ and B Xs£~^£~ . The lead- 
ing non-perturbative corrections to the free quark decay model, of relative order 1/m^, may 
be expressed entirely in terms of two low-energy parameters. One of these is determined 
from the splitting between the heavy pseudoscalar and vector mesons; a model-dependent 
estimate of the other comes from QCD sum rules. In addition to the total rates, we have 
computed the correction to the average photon energy in i? ^ Xs'j and found the shift to 
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be small. The correction to the spectrum for B Xsl~^l~ is larger and for high invariant 
mass lepton pairs is at about the 10% level. 

Finally, we note that there has been considerable recent work in which a similar 
formalism has been applied to semileptonic b decays . 
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Appendix A. Matrix Elements and the Equation of Motion 

In this appendix we derive an identity for the matrix element of a time-ordered product 
of two operators, where one of the operators vanishes by the equation of motion of the 
theory. This will be a generalisation of a proof by Politzer that matrix elements of 
single operators which vanish by the equation of motion themselves vanish. We will derive 
our result within the context of the heavy quark effective theory (HQET), because this 
is the application which we have in mind, but with obvious modifications our result is 
completely general. 

We begin by recalling how such time-ordered products arise within HQET. In this 
effective theory, the heavy quark part of the lagrangian takes the form 

C = hviDh + ^hF{D)h + ..., (A.l) 
zrub 

where 

F(D) = (iD)2 - Zi{p){v ■ iD)2 + Z2{^i)s^''G^, , (A.2) 

and the ellipses denote terms of higher order in the l/m^ expansion. Here the gluon field 
strength is defined by G^i, = [iD^, iD,^], and s^'^ = —^a^^. The renormalisation constants 
Zi[n) and Z2{ij.) are given in Section 2. The equation of motion in HQET is derived from 
the leading term in the lagrangian ([A.l| ), and is simply 



v-iDh = 0. 
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(A.3) 



Instead of being included in the equation of motion, the corrections to L in eq. ( |A.1| ) are 
treated as perturbations. They reappear in the following way: because the states \M) 
of HQET are defined by the truncated equation of motion ( |A.3| ), they differ from those 
\B) of full QCD. The states |M) have the significant advantage that, unlike they are 
independent of the heavy quark mass m^, and so have simple transformations under the 
spin- flavor symmetries of the effective theory. The difference between |M) and the physical 
states \B) is then compensated by including in the matrix elements of effective operators 
additional time-ordered products with the subleading terms in L That is, if we 



have an effective operator h C(D)h whose matrix element we require between eigenstates 
l-B) of full QCD, then we must write 

{B{p')\h'c{B)h\B{p)) = {M{v')\h'CiB)h\M{v)) 

+ ^ i [ dx {M{v') I T{h'C{D)h, hF{D)h{x) }\M{v)) + ... . 

(A.4) 

We have shown the expansion up to order l/rrib explicity; the ellipses denote denote terms 
of higher order which may be included if more accuracy is needed. We consider here the 
general case in which the initial and final heavy quarks have different four-velocities. The 
field h creates a heavy quark with velocity v^, while h' creates one with velocity v'^. There 
is a separate effective lagrangian ( |A.1| ) for each of these fields, but for simplicity we will 
include the l/rrib corrections only for the field h. The time-ordered products in ( [A.4| ) are 
new nonperturbative matrix elements which must be evaluated if one wishes to use the 
effective theory beyond leading order. 

We will be concerned with a special case of eq. ( |A.4[ ), in which the operator h C{D)h 



vanishes by the equation of motion ( [A.3| ) of the effective theory.0 That is, C(D) takes the 
particular form 

C(D) = A(D)^;-iD, (A.5) 

where A(D) may include an arbitrary Dirac structure. Politzer has shown that matrix 
elements of such operators, such as would appear in the first term in eq. ( |A.4| ), vanish 

{M{v')\h'A(D)v ■ iDh \M{v)) = . (A.6) 



^ We are grateful to A. Manohar for discussions of this point. See also ref. [20| 
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Note that it is the effective theory states \M) which appear here. The purpose of this 
appendix is to generahse this argument to prove a similar identity for the time-ordered 
product appearing in the second term of ( |A.4| ) , namely, that 



i j dx{M{v')\T{h! A{I^)v ■iDh,hF{T>)h{x)]\M{v)) ^^^^ 

= - {M{v')\t A{V))P+F{V))h\M{v)) , 

where P-f = i(l + ^). For the computation of this paper, we will apply this identity in 
the case that A(D) actually contains no derivatives, and at zero recoil; then the matrix 
element on the right-hand side of (|A.7D is of the simple form 



(M|/iriD^iD^/i|M), (A.8) 

and can be evaluated in terms of the constants Ai and A2 as in Section 2. 

Before proving the identity ( [A.7|) , however, we note its relation to the result we would 
obtain by taking a different approach. Instead of introducing the states |M), which are 
eigenstates of the lowest order effective lagrangian, we could chose to work always in terms 
of the full states \B) of QCD. This would be undesireable, in that it would reintroduce 
the mass-dependence which it is the goal of HQET to remove, thereby obscuring the spin- 
flavor symmetries of the heavy quark limit. However, if we do so, the equation of motion 
is given by the full Lagrangian ( [A.l| ), taking the form 

viT)h = 3— P+F(D)/i + . . . . (A.9) 

This equation of motion may be applied directly to matrix elements between the states 
l-B). We then flnd for the matrix element ( |A.4]) the relation 

{B{p')\h!c{Ji)h\B{p)) = --^{B{p')\h'AiD)P+F{D)h\B{p)). (A.IO) 

Zrub 



Inserting eqs. ( |A.6|) and (|A.7| ) into eq. ( |A.4] ), and noting that the states \M) and \B) 



differ only at order l/rrii,, we see that this is the same result that we flnd working entirely 
within the effective theory. The proof which we now present may be seen as verifying the 
consistency of the effective theory approach. It is both an application and a generalisation 
of the proof of Politzer. 
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We begin by using the LSZ reduction formula to write the desired matrix element in 
terms of a vacuum expectation value: 

i J dx {M{v')\ T{h'A{'D)v ■ iDh{0), hF(D)h{x)}\M{v)) 

= ijdxj dzdz'e'^-'e'^'-''v ■ iD^ v' ■ (^-H) 

• (0| T{qr'h\z'), h!A{T>)v ■ iD/i(0), hF{T>)h{x), hVq{z)]\Q) . 

Here the operators /iFq and qF'/i interpolate the initial and final meson states, respectively. 
(Of course, the proof is valid for any external heavy hadrons, not just mesons.) Note that 
in HQET, the one-particle poles are projected out by the differential operator v ■ iD rather 
than by [(iD)^ - M|] as in fuU QCD. 

We now write the generating functional for Green's functions of this theory: 

exp(iW^) = j [dh] [dh] [dh'] [dt] [d^^] exp |i y dy [£o + + + -Sm + . . .] I . (A.12) 

Here 

Co^hv- iDh + h'v' ■ iDh' (A. 13) 

is the lagrangian of the effective theory, and we have included explicitly a variety of relevant 
source terms: 

Sj = Jh'A{D)viDh, 

SL = LhF{D)h, (A.14) 
Sm = KhVq + K'qT'h' . 

The ellipses denote sources for the fermions and gauge fields, and gauge-fixing terms which 
will play no role in the analysis. With these definitions, then, we have 

exp(iVy ) 



SiK'iz') 5iJ{0) 5iL{x) 5iK{z) ^ ^ ' ,,„rces=o (A.15) 
= (0| T{qr'/i'(^')> ^A{p)v ■ iD/i(0), }iF{p)h{x), /^rq(^)}|0) . 

We now perform a shift of the integration variable /i, 

h = TC - jtA{p)P+, (A.16) 
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insert it into the generating functional (|A.12|) , and drop terms of order . We then obtain 
shifts in some of the expressions ( |A.13| ) and ( |A.14|) : 

Cq = h*v ■ iDh - Jh'A{D)v ■ iDh + h'v' ■ iDh' , 
Sj = Jh'A(D)viDh, 

(A.17) 

SL = Lh F{D)h - LJ h A(D)P+F{D)h , 
Sm = Kh*Tq - KJh'A{D)P+Tq + K'qT'h' . 

Note that the original source term Sj cancels against the shift in Co in (|A.17| ), but new 
terms appear in 5"^ and Sm- Replacing the dummy variable h by h, we recover the 
generating functional ( |A.12D , but with the source Sj changed, 



Sj -LJ h A{'D)P+F(D)h - KJ h A{D)P+rq. (A.18) 



We now repeat the derivative in eq. (|A.15| ). When we set the sources to zero, we see that 



a derivative with respect to L or K must come with a derivative with respect to J to give 
a nonzero contribution. We then find 

{0\T{qr'h\z'),h'A{D)viDh{0), hF{D)h{x), 7^rq(z)}|0) 

= i5{z) (0| T{qT'h'{z'), hA{D)h{x), h' A{D)P+Tq{0)}\0) (A.19) 
+ i6{x) (0| T{qr'h\z'), h'A{D)P+F{D)h{0), hrq{z) } |0) . 

Finally, we must perform the integral {i f dx f dzdz'e^^'^e^^ v ■ iD^f' ■ iDz') to re- 
cover the matrix element ( A.ll|) . In this integral, the first term on the right-hand side of 



eq. ( |A.19| ) vanishes for an on-shell state with v ■ k = 0, because the integral over z is trivial 
and there is no longer a one-particle pole to pick out. The second term, however, yields 
an iS'-matrix element in the usual way, 

i/d.i.w/d.d.'e-e--..iD.„'.iD. 

■ {0\T{qT'h'{z'), h'A{D)P+F{D)h{0), hTq{z)}\0) (^-20) 

= -{M{v')\h' A{D)P+F{D)h{0) \M{v)) . 

Thus we obtain the desired identity, 

i [ dx{M{v')\T{h'A{D)viDh,hF(D)h{x)]\M{v)) 
J (A.21) 

= - {M{v')\h'A{D)P+F{'D)h\M(v)) . 
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Note that the term in -F(D) proportional to {v ■ iD)^ will not contribute here, since this 
matrix element is of the form (|A.6| ) and hence vanishes by the equation of motion. 

A few additional comments are in order. First, Politzer's result ( |A.6|) for the matrix 
elements of an operator which vanishes by the equation of motion follows from an identical 
derivation, but with the derivative d/6iL omitted. In this case the second term of eq. ( |A.19| ) 
does not appear, and we obtain zero instead of the right-hand side of eq. ( |A.2U| ) . (We stress 
that the intermediate result ( |A.19| ), which is the key to both proofs, is derived by Politzer 
in full generality.) Second, it is clear how this result is to be generalised to the time-ordered 
product of an arbitrary number of operators. Essentially, we obtain a term on the right- 
hand side for each contraction of h A(D)v ■ iDh with an operator insertion hG{D)h, where 
G(D) is any function of covariant derivatives. More than one contraction may be required 
to give a nonzero result; for example, an operator of the form h A(D){v ■ iD)'^h will have a 
nonvanishing matrix element only when included in a time-ordered product with n other 
operators such as hF(D)h. Finally, we reiterate that while for concreteness we have framed 
our derivation within the heavy quark effective theory, it is in fact completely general. 
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Figure Captions 



Fig. 1. Feynman diagrams contributing to the time-ordered product Tjbris, sr2b}. 

Fig. 2. Invariant mass spectrum for B — > Xg£'^£~ . The sohd hne corresponds to the 
parton model, while the dashed line corresponds to Ai = 0.5 GeV^ and A2 = 
0.12 GeV^. The cusp at s = {2mc/mhf' corresponds to the charm threshold. 
Near this point our estimate of the nonperturbative corrections is not valid due 
to resonance effects. 
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